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AN ENERGY-BASED HYSTERESIS MODEL FOR MAGNETO STRICTIVE 

TRANSDUCERS 

F.T. CALKINS*, R.C. SMITH*, AND A.B. FLATAU* 

Abstract. This paper addresses the modeling of hysteresis in magnetostrictive transducers. This is 
considered in the context of control applications which require an accurate characterization of the relation 
between input currents and strains output by the transducer. This relation typically exhibits significant 
nonlinearities and hysteresis due to inherent properties of magnetostrictive materials. The characterization 
considered here is based upon the Jiles- Atherton mean field model for ferromagnetic hysteresis in combination 
with a quadratic moment rotation model for magnetostriction. As demonstrated through comparison with 
experimental data, the magnetization model very adequately quantifies both major and minor loops under 
various operating conditions. The combined model can then be used to accurately characterize output strains 
at moderate drive levels. The advantages to this model lie in the small number (six) of required parameters 
and the flexibility it exhibits in a variety of operating conditions. 

Key words. Hysteresis model, magnetostrictive materials 

Subject classification. Applied and Numerical Mathematics 

1. Introduction. This paper addresses the modeling of hysteresis in magnetostrictive transducers. 
The capabilities for actuation and sensing in such transducers are provided by the dual magnetostrictive 
effects in the core material: (i) the application of a magnetic field generates strains in the material and (ii) 
material stresses yield measurable magnetic effects. One core magnetostrictive material which has proven 
very effective at room temperatures and nominal operating conditions is Terfenol-D (see [1, 2] for descriptions 
of the material and its capabilities). Due to the magnitude of the strains and forces generated by the material, 
Terfenol-D transducers have been employed as ultrasonic transducers, sonar projectors and provide the 
capability for controlling vibrations in heavy structures and industrial machinery. 

Several properties inherent to magnetostrictive materials must be addressed when designing systems 
which employ them. The first concerns the hysteresis and nonlinear dynamics exhibited by the materials. 
This is due to inherent magnetic properties of the materials and is particularly pronounced at higher drive 
levels. It is also well documented that Terfenol-D performance is highly sensitive to operating conditions 
such as temperature, mechanical prestress, magnetic excitation (bias and AC amplitude), frequency and 
external load [3, 4, 5]. Several of these aspects (e.g., prestress and external loads) involve system aspects 
external to the core Terfenol-D material which makes the extrapolation of results from isolated laboratory 
samples to actual transducer design difficult and motivates consideration of the transducer as a whole. 

Accurate modeling of transducer dynamics is necessary to take advantage of the full capabilities of 
the materials and to provide the ability for tailoring the performance of the transducers by modifying easily 
adjusted operating conditions. To attain these objectives, the model must accurately characterize both major 
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(symmetric) and minor (nested and asymmetric) hysteresis loops as well as constitutive nonlinearities. The 
model must also incorporate the sensitivities with respect to operating conditions and be in a form amenable 
for eventual incorporation in models for underlying structural systems. Finally, the model must be suitable 
for controller design in the sense that it is efficient to implement and characterizes all dynamics which may 
be specified by the control law. For example, a model which characterizes major loops but not minor ones 
would be less useful in a feedback control law which cannot differentiate between the two. 

The model we consider is obtained through the extension of the ferromagnetic mean field model of Jiles 
and Atherton [6, 7, 8, 9] to magnetostrictive transducers. This provides a characterization for the inherent 
hysteresis which is based upon the anhysteretic magnetization along with reversible and irreversible domain 
wall movements in the material. When coupled with nonlinear strain/magnetization relations, this yields a 
model which characterizes strain outputs in terms of input currents to the driving solenoid. Minor loops arc 
incorporated through the enforcement of closure conditions. 

With regard to design criteria, this model is currently constructed for a transducer with quasi-static input 
and fixed temperatures (these are commonly employed conditions for initial transducer characterization). 
The capability for having different prestresses and variable input magnitudes to the driving solenoid are 
included in the model and demonstrated through comparison with experimental data. The advantages 
of this approach lie in the accurate fits attainable in the considered regimes with a small number (six) 
of physical parameters to be identified through least squares techniques. This provides the method with 
significant flexibility and low computational overhead. The model is also in a form which can be extended to 
variable temperature and frequency regimes and can be incorporated in a large variety of structural models 
(e.g., see [10, 11]). As a result, it shows great promise for use in transducer design for precision positioning 
and structural controllers [12]. 

To place this modeling approach in perspective, it is useful to briefly summarize existing techniques 
for characterizing magnetostrictive transducers. For initial applications, linear field/magnctization relations 
were used to approximate the transducer dynamics [1, 13]. While this approach is reasonable at low drive 
levels, it is inaccurate at moderate to high input levels due to inherent hysteresis and material nonlinearities. 
In this latter regime, various phenomenological or empirical techniques, including Preisach models, have been 
employed to quantify the input/output relations [14, 15]. Phenomenological approaches circumvent unmod- 
eled or unknown physical mechanisms and have the advantage of generality. While some connections have 
been made between underlying physical processes and Preisach models [16], this genre of model typically pro- 
vides less insight into physical dynamics than a model developed from physical principles. Furthermore, such 
empirical models generally require a large number of nonphysical parameters and are not easily adapted to 
changing operating conditions. This increases implementation time [17] and will limit flexibility if employed 
in a control law. 

A typical magnetostrictive transducer is described in Section 2. This illustrates the system being modeled 
and indicates design issues which must be incorporated in the model. The energy-based model is discussed 
in Section 3 and the applicability of the model in a variety of experimental settings is presented in Section 4. 
These results illustrate the accuracy and flexibility of the model at fixed temperatures and low frequencies 
and indicate the extensions necessary for use in other regimes. 

2. Magnetostrictive Transducers. The issues which must be addressed when developing a compre- 
hensive model are illustrated in the context of the transducer depicted in Figure 1. As detailed in [14], this 
construction is typical for actuators currently employed in many structural applications; hence it provides 
a template for the development of models which will ultimately enhance design and performance. Details 
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regarding the specific experimental setup used here are provided in Section 4. 

From a design perspective, the transducer can be considered as the entire system which facilitates the 
utilization of the magnetostrictive core for applications. For modeling purposes, the key components arc 
the magnetostrictive core, a DC magnetic circuit, a driving AC circuit and a prestress mechanism. The 
magnetostrictive material used in the transducer for the experiments reported in Section 4 was comprised 
of Terfenol-D, Tbo.sDyo.yFei.g, while the driving AC magnetic field was generated by a surrounding wound 
wire solenoid. As illustrated by the experimental data plotted in Figure 2, the relationship between the 
applied field H and resulting magnetization M exhibits significant hysteresis while the relationship between 
the magnetization and strain e is highly nonlinear. Moreover, the strains in an unbiased rod are always 
positive since the rotation of moments in response to an applied field always produce an increase in length. 
To attain bidirectional strains, a DC bias is provided by the enclosing cylindrical magnet (alternatively, a 
biasing DC current could be applied to the solenoid). Finally, the prestress bolt further aligns the orientation 
of magnetic moments and maintains the rod in a constant state of compression. 

To fully utilize the transducer for structural applications and eventual controller design, it is necessary 
to characterize the relationship between the current I applied to the solenoid, the resulting field H, the 
associated magnetization M and finally, the generated strains e. A characterization based upon the Jiles- 
Atherton ferromagnetic hysteresis model is presented in the next section. 



Wound Wire Solenoid 

Figure 1 . Cross section of a typical Terfenol-D magnetostrictive transducer. 



Magnetic Reid (H) x 10 4 Magnetization (M) x to 5 


(a) (b) 

Figure 2. Relationship in experimental data between (a) the magnetic field H and the magnetization A/, 
and (b) the magnetization M and the generated strains e. 
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3. Domain Wall Dynamics. The transducer model described here is based upon the theory that 
magnetization in ferromagnetic materials is due to the realignment of magnetic moments within the material. 
Such materials exhibit the property that at temperatures below the Curie point , moments are highly aligned 
in regions termed domains (the reader is referred to [6, 18] for further discussion regarding the experimental 
verification of domain properties) . The reorientation of moments can occur both in bulk within the domains 
or within transition regions, termed domain walls, between domains. 

For a material which is defect free, the former mechanism leads to anhysteretic (hysteresis free) be- 
havior which is conservative and hence reversible. Such a situation is idealized, however, since defects are 
unavoidable (e.g., carbides in steel) and in many cases, incorporated in the material to attain the desired 
stoichiometry (e.g., second-phase materials such as Dysprosium in Terfenol-D). These defects or inclusions 
provide pinning sites for the domain walls due to the reduction in energy which occurs when the domain 
wall intersects the site. For low magnetic field variations about some equilibrium value, the walls remain 
pinned and the magnetization is reversible. This motion becomes irreversible at higher field levels due to 
wall intersections with remote inclusions or pinning sites. Note that pinning effects lead to phenomena such 
as the Barkhausen discontinuities observed in experimental magnetization data [6, 18]. The energy loss due 
to transition across pinning sites also provides the main mechanism for hysteresis in ferromagnetic materials. 


Magnetostriction 


The model presented here ultimately provides a relationship between the current I input to the solenoid 
and the strain e output by the transducer. As a first step, we characterize the magnetostriction which 
results at a given magnetization level. The magnetostriction A = ^ indicates the relative change in length 
of the material from the ordered, but unaligned state, to the state in which domains are aligned. While the 
magnetostriction does not quantify DC effects, the effects of domain order, or thermal effects, it does provide 
a measure of the strains generated in a Terfenol transducer. 

As detailed in [6] , consideration of the potential energy for the system yields 


(3.1) 


_ 3 A s 2 
2 M? M 


as a first approximation to the relationship between the magnetization and magnetostriction. Here M s 
and X s respectively denote the saturation magnetization and saturation magnetostriction. For an isolated 
Terfenol-D sample, M s represents the magnetization required to rotate all moments and has been observed 
to have the approximate value M s ~ 7.9 x 10 5 A/m [19]. This parameter has a similar interpretation in 
the full transducer model but will be shown in the examples of the next section to have the slightly smaller 
value of M s = 7.65 x 10 5 A/m. This illustrates the necessity of estimating such parameters for the specific 
transducer under consideration. The value of A* depends upon the initial orientation of moments and hence 
upon the applied prestress. In the absence of applied stresses and under the assumption of a cubic anisotropy 
model, Ag can be defined in terms of the independent saturation magnetostrictions Aioo and Am in the (100) 
and (111) directions, respectively. As detailed in [6], under the assumption that the material contains a large 
number of domains and has no preferred grain orientation, averaging of domain effects yields the expression 


X s = -Aioo + 5^111 

for the total saturation magnetostriction (typical saturation values for Terfenol arc Aioo = 90 x 10 -6 and 
A nl = 1600 x 10 -6 ). As will be noted in the examples of the next section, this saturation value is highly 
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dependent upon the operating conditions (e.g., applied prestress) and the parameter X s must be estimated 
through least squares techniques for the specific conditions under consideration. 

For the operating conditions under consideration, the quadratic expression (3.1) adequately models the 
relationship between the magnetization and strain at low to moderate drive levels. For higher drive levels and 
frequencies along with variable temperature and stress conditions, however, it must be extended to include 
higher-order mechanisms (e.g., the sensitivity of the system with regard to changing stress is an important 
and well documented phenomenon [3, 4, 5]). This can be accomplished through the incorporation of stress 
dependence in A s and the use of higher-order magnetostrictive models as discussed in [20]. Alternatively, 
higher-order effects and magnetostrictive hysteresis can be incorporated through an energy formulation 
as detailed in [9]. Finally, the effects of magnetomechanical coupling and mechanical resonances must be 
incorporated in various operating regimes. Hence this component of the transducer model should be extended 
as dictated by operating conditions. 

We next turn to the characterization of the magnetization M in terms of the input current /. To 
accomplish this, it is necessary to quantify the effective field H e fj associated with the magnetic moments in 
the core material, the anhysterctic magnetization M an , the reversible magnetization M rev and the irreversible 
magnetization M{ rr . 

Effective Magnetic Field 

In general, the effective field is dependent upon the magnetic field generated by the solenoid, magnetic 
domain interactions, crystal and stress anisotropies, and temperature. For this model, we are considering the 
case of fixed temperature and compressive prestresses in excess of 0.8 ksi. It is noted in the computations on 
pages 126 and 410 of [6] that for polychrystaline Terfenol, a compressive stress of a = 6.25 MPa or 899 psi 
is required to align moments perpendicular to the stress under the assumptions that X s — 1067 x 10 ~ 6 and 
the crystal anisotropy constant is K\ — — 2 x 10 4 J/m 3 . While this computed value of <7 is highly dependent 
upon temperature and operating conditions, it indicates that stress anisotropies will start to dominate crystal 
anisotropies by 1 ksi with the effect magnified at higher prestresses. This motivates the use of a model which 
neglects crystal anisotropies when operating in high stress regimes. 

Under the assumption of fixed temperature and sufficiently large prestresses, the effective magnetic field 
is modeled by 


H eff = H + aM + H a 


where H — nl is the field generated by a solenoid with n turns per unit length, aM quantifies the field due 
to magnetic interactions between moments, and H a is the field due to magnetoelastic domain interactions. 
The parameter a quantifies the amount of domain interaction and must be identified for a given system. 
The field component due to the applied stresses can be quantified through thermodynamic laws to obtain 


H a 


3 a 
2 Ho 



a ,T 


(see [9, 20] for details). Here fi o is the free space permeability, and the subscript T denotes constant 
temperature in degrees Kelvin. Note that with the approximation (3.1) for A, the effective field can be 
expressed as 


H ef j = H + aM. 


where a = a + § -r*£r 


Mo 
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Anhysteretic Magnetization 

The anhysteretic magnetization is computed through consideration of the thermodynamic properties of 
the magnetostrictive material. Under the assumption of constant domain density AT, Boltzmann statistics 
can be employed to yield the expression 

(3.2) M an = M s C(H e ff/a) 

where C(z) = coth(z) - l/z is the Langevin function. The constant a is given by a = where k B 

is Boltzmann’s constant and k B T represents the Boltzmann thermal energy. We point out that a cannot 
directly be computed for a transducer due to the fact that N is unknown. Hence it is treated as a parameter 
to be identified for the system. We also note that this expression for M an is valid only for operating conditions 
under which H e jf is valid. For example, if prestresses are sufficiently small so that crystal anisotropies are 
significant, the expression must be modified to incorporate the differing anisotropy energies in the different 
directions. One approach to modeling the effects of anisotropy is given in [21], 

Irreversible , Reversible and Total Magnetization 

The anhysteretic magnetization incorporates the effects of moment rotation within domains but does not 
account for domain wall dynamics. As noted previously, the consideration of domain wall energy yields addi- 
tional reversible and irreversible components to the magnetization. The consideration of energy dissipation 
due to pinning and unpinning of domain walls at inclusions yields the expression 

dMi rr M a n M{ rr 

( ' J ~dH~ ~k6^5 (M„ n - M irr ) ^ 

for the differential susceptibility of the irreversible magnetization curve [7, 20]. The constant k = 2 ^}(i-c) ' 
where (p) is the average density of pinning sites, (e n ) is the average energy for 180° walls, c is a reversibility 
coefficient, and m is the magnetic moment of a typical domain, provides a measure for the average energy 
required to break a pinning site. The parameter S is defined to have the value +1 when > 0 and —1 
when < 0 to guarantee that pinning always opposes changes in magnetization. In applications, 5 can 
be directly determined from the magnetic field data while k is identified for the specific transducer and 
operating conditions. 

The reversible magnetization quantifies the degree to which domain walls bulge before attaining the 
energy necessary to break the pinning sites. As derived in [7], to first approximation, the reversible magne- 
tization is given by 

(3.4) M rev — c(M an Mirr') • 

The reversibility coefficient c can be estimated from the ratio of the initial and anhysteretic differential sus- 
ceptibilities [8] or through a least squares fit to data. Properties of all the model parameters are summarized 
in Table 1. 

The total magnetization is then given by 

(3.5) M = M rev + Mi rr 

with M irr and M rev defined by (3.3) and (3.4) and the anhysteretic magnetization given by (3.2). The full 
time-dependent model leading from input currents to output magnetization is summarized in Algorithm 1. 
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When combined with (3.1), this provides a characterization of the output strains in terms of the current I in- 
put to the solenoid. Note that this model is valid for fixed temperature and quasi-static operating conditions. 
The extension to more general operating conditions will involve the previously mentioned modifications to 
the effective field. 


(i) H(t) = nl(t) 


(ii) H eff (t) = H(t) + aM(t) + H a {t) 


(iii) 


M an {t) - M s 


coth 


V a ) \H eff (t) 


\ dM irr , v ^ dd M an {t) Mi Tr {t^ 

~ n kS - 5[M on (t) - M, rr (f)]^ 

(v) M rev {t) = c[M an (t) - M irr {t)} 


(vi) M(t) = M rev (t ) + M irr (t) 


Algorithm 1 . Time- dependent model quantifying the output magnetization M(t ) in terms of the input 
current I(t). The parameter a is given by a = a + | where <7 q is the applied prestress. 


Parameter 

Physical Property 

Effects on Model 

a 

Quantifies domain in- 
teractions 

Increased values lead to steeper slopes for an- 

hysteretic and magnetization 

curves. 

a 

Shape parameter 
for M an 

Increased value decreases slope of M an . 

k 

Average energy 
required to break pin- 
ning sites 

Increased value produces wider hysteresis 
curve and narrower minor loop. 

c 

Reversibility coef- 
ficient 

Decrease in value leads to wider hysteresis 
curve. 

M s 

Saturation 

magnetization 

Increase leads to large saturation value for 
magnetization. 

K 

Saturation 

magnetostriction 

Increase leads to large saturation value for 
magnetostriction. 


Table 1 . Physical properties and effects of model parameters a,a,k,c,M s , X s . The parameter a is then 
given by a = a + | where <7 0 is the apphed prestress. 
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Asymmetric Minor Loops 


The final aspect which we consider here concerns the modification of the model to incorporate minor 
(asymmetric) loops. Such loops occur when the sign of ^ is reversed for a trajectory lying within the 
interior of the major loop. To preserve order in the sense that forward paths do not intersect, it is necessary 
that minor loops close. The model (3.5) can be employed for the first half of the minor loop but does not 
ensure closure. This property is incorporated in the model through the consideration of a working volume 
and volume fraction for either the magnetization or the reversible and irreversible components. 

To illustrate the first case, we let to, ti and t 2 respectively denote the times when the minor loop starts, 
when it turns due to a change in the sign of and when it closes (see Figure 3). The corresponding values 
of the magnetic field and magnetization are H(to), H (£i), H(t 2 ) and M(£o), M(t 2 ). Note that in order 

to guarantee closure of the minor loop, it is necessary to require that H(to) = H(t 2 ) and M(t 0 ) = M(t 2 ). 
Direct integration of (3.5) yields 

J f t2 dM 

~ T~ ds 

ti ds 


which in general will not be equal to M(to). To attain closure, we define 


M{t) = M(ti) - 


M(ti)-M(* 0 ) [ l dM 

M 2 -M(t l ) J tl 


ds 


-ds 


for t e [t\ } t 2 \. The magnetization values M(t 0 ), M 2 and are computed using (3.5). Through the 

inclusion of this volume fraction 


(3.6) 




the magnetization is forced to satisfy the closure property M(t 2 ) — M(£o)- A similar formulation of volume 
fractions for the component reversible and irreversible magnetizations is given in [22] while extensions of the 
model to accommodate more complex anhysteretic effects can be found in [23]. 

The viability of the model with minor loops closed via (3.6) is illustrated in the next section. We note 
that for the operating conditions targeted in this paper, the model accurately characterizes the transducer 
response including both major loops and nested minor loops. 


/ (H(t 2 ), M 2 ) 



Figure 3. Closure requirements for minor loops. 
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4. Model Fits to Experimental Data. The model fits to experimental transducer data using the 
relations summarized in Section 3 are presented here. Following a description of the experimental transducer, 
two cases are considered. The first illustrates the performance of the magnetization and magnetostriction 
models under various drive levels with a 1.3 ksi prestress applied to the Terfenol rod. Included in these 
results are model fits to data which contains minor loops. The second case illustrates the performance of 
the model for a prestress of 1.0 ksi. As discussed in the last section, the stress-dominated anisotropy model 
for the magnetization is valid for both cases. Taken in concert, these examples illustrate the accuracy and 
flexibility of the magnetization model for a range of drive levels, magnetic biases and prestresses for quasi- 
static operating conditions at fixed temperature. The quadratic magnetostriction model is also accurate at 
low to moderate drive levels but must be extended to incorporate the hysteresis and saturation present at 
high drive regimes. 

Experimental Transducer 

The experimental data reported here was collected from a broadband Tcrfcnol-D transducer developed 
at Iowa State University. The nominal resonance range was designed for structural applications (1-10 kHz). 
Furthermore, the transducer was designed to produce an output free from spurious resonances and to permit 
adjustable prestress and magnetic bias. 

The Terfenol-D (Tbo.3Dyo.7Fe1. 9) rod employed in the transducer had a length of 115 mm and a 12.7 mm 
diameter. The rod was placed inside two coils consisting of an inner single layer 110-turn pickup coil and a 
multi-layer 800-turn drive coil. A current control amplifier (Techron 7780) provided the input to the drive 
coil to produce an applied AC magnetic field and DC bias as necessary. The reference signal to this amplifier 
was provided by a Tektronix spectrum analyzer and the applied magnetic field H generated by the drive coil 
had a frequency of 0.7 Hz and magnitude up to 5.6 kA/m (700 Oe ) per ampere. The pickup coil was used 
to measure the induced voltage from which the time rate change of the magnetic induction B was computed 
using the Faraday- Lenz law. 

A cylindrical permanent magnet surrounding the coils provided the capability for generating additional 
DC bias if necessary. This permanent magnet was constructed of Alnico V and was slit to reduce eddy current 
losses. Note that for the experiments reported here, biases generated in this manner were unnecessary and 
the reported data is unbiased (i.e. , the permanent magnet was demagnetized). Finally, mechanical prestresses 
to the rod were generated by a variable prestress bolt at one end of the transducer and Belleville washers 
fitted at the opposite end of the rod. 

The measurable output from the transducer included the current and voltage in the drive coil, the 
voltage induced in the pickup coil, and the mechanical output. To quantify the mechanical output, a Lucas 
LVM-10 linear variable differential transformer based upon changing reluctance was used to measure the 
displacement of the transformer output interface connection. Temperature was maintained within 5° C of 
the ambient temperature (23° C) by monitoring two thermocouples attached to the Tcrfenol-D sample. 

Parameter Estimation 

The use of the magnetization and magnetostriction models to characterize transducer dynamics requires 
the estimation of the parameters a, a, fc,c, M s and \ s summarized in Table 1. The parameters 5, a, A: and 
c are in essence averages which arise when extending physics at a microscopic level to the macroscopic 
scale necessary for control implementation. Hence, while they have physical interpretations and tendencies, 
they must be estimated for individual transducers. The parameters M s and \ s are macroscopic and have 
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published values for Terfenol under various operating conditions. Sufficient variation occurs in the values, 
however, that we also estimated them for the individual transducer. 

The full set of parameters was estimated through a least square fit with experimental data from the 
previously described transducer. The optimization was performed in two steps. In the first, the values of 
q = (5,0,0, k , M s ) were estimated through minimization of the functional 

s 

(4.1) •/(*) = £ \M(t i ;q)-z i \ 2 

i—1 

where z x denotes the experimentally measured value of the Terfenol magnetization at time t x . The modeled 
magnetization at time t t for parameter values q is denoted by M(ti ; q) (see (3.5) or (vi) of Algorithm 1). The 
functional (4.1) was minimized using a constrained optimization algorithm based upon sequential quadratic 
programming (SQP) updates. 

With the estimated values of 5, a, c, k and M a , the model fits to the experimental magnetization curves 
can be obtained. The second step concerns the estimation of X s to attain reasonable fits in the magnetostric- 
tion model (3.1). This was accomplished through a least squares fit with displacement or strain data from 
the transducer. 

Initial magnetization parameters were estimated using this technique for the transducer with an applied 
prestress of 1.3 ksi. The resulting values are summarized in Table 2 while model fits are illustrated in 
Figure 4. From strain data, the saturation magnetostriction constant for this case was determined to be 
X s = 1003 x 10 -6 for high drive levels and A s = 1221 x 10” 6 at low drive levels (the difference in values is 
further discussed in the next section). 

To ascertain the robustness of the model with respect to applied prestresses, we then considered the 
estimation of parameters and performance of the model with a prestress of 1.0 ksi . For this case, we 
fixed the parameters a, c, A I s which have the least theoretical dependence upon prestress and estimated the 
parameters A:, 5, X a through a least squares fit to the data. The estimated magnetization parameters are 
again summarized in Table 2 while the saturation magnetostriction was found to be X s = 995 x 10 ~ 6 at high 
drive levels. 

A comparison of the estimated values of k indicates significant changes due to the effects of stress on 
the pinning energy at magnetic inclusions. The change in the saturation magnetostriction X s is due to 
stress-induced changes in the initial domain configuration. The stress dependence in 5 = a + § is 
primarily due to magnetomechanical stress anisotropies which are quantified by the term | • Note that 

for the compressive prestress <7o = —1 ksi ~ —6.9 MPa and the estimated values for A a , a, M sy the magnetic 
coupling parameter a has the computed value a = 0.032 while it has the value a = 0.035 for <7o — —1.3 ksi . 
This small variation in the values of a (less than 9%) illustrates the consistency of the model with regard to 
nearly constant applied stresses. Moreover, it indicates that one has the capability for identifying and fixing 



(To = 1-3 ksi 

00 = 1.0 ksi 

M s {A/m) 

7.65 x 10 5 

7.65 x 10 5 

a (A/m) 

7012 

7012 

c 

0.18 

0.18 

k (A/m) 

3942 

3283 

a 

-0.02 

-0.01 


Table 2. Estimated magnetization parameters for the transducer with prestresses of 1.3 ksi and 1.0 ksi. 
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the parameter a and incorporating subsequent stress effects through the component | • The use of this 

strategy has been substantiated by the highly accurate model fits obtained with fixed a. 

Magnetization Model 

We consider first the performance of the quasi- static magnetization model summarized in Algorithm 1 
under a variety of operating conditions. The model is formulated to be flexible with regard to various drive 
levels and prestresses and it was within this regime that the performance was tested. Data was collected 
at multiple drive levels with prestresses of 1.0 ksi and 1.3 ksi applied to the Terfenol rod. As detailed in 
[3], prestresses within this range yield nearly optimal magnetomechanical coupling and strain coefficients 
for the specific transducer. Parameters for the magnetization model were estimated through the previously 
described least squares techniques and used to obtain model responses under the various conditions. In each 
case, the measured applied field H was used as input to the model. 



-« -4 -2 0 2 4 6 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 



(c) (d) 


Figure 4. Experimental data ( ) and magnetization model dynamics ( ) for multiple drive and 

prestress levels; (a) Three drive levels with 1.0 ksi applied stress, (b) Magnified view of 1.0 ksi case, (c) Two 
drive levels with 1.3 ksi applied stress, (d) Magnified view of 1.3 ksi case. 
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The model fits at three drive levels for the 1.0 ksi case arc illustrated in Figure 4a, b while fits for 
two drive levels with a 1.3 ksi applied stress are illustrated in Figure 4c, d. For each fixed prestress, the 
same fixed parameters in Table 2 were used to attain the model responses at the multiple drive levels. The 
variation in model dynamics is due solely to the changes in the input fields. This illustrates the flexibility 
of the model with respect to drive levels. As noted in previous discussion and summarized in Table 2, only 
the parameter k and stress contribution § a mus t be modified to account for changes in prestress. 

Hence the model is also highly flexible with respect to applied prestresses. 

Close examination of Figure 4a, c indicates that one aspect of the experimental transducer behavior 
which is not quantified by the model is the constricted or ‘wasp-waisted’ behavior which occurs at low 
applied fields. This behavior has been noted by other investigators [18, 24] and is hypothesized to be due 
to 180° domain changes [25]. While quantification of this effect is ultimately desired, the accuracy and 
flexibility of the current magnetization model are sufficient for control applications in this operating regime. 

Magnetostriction Model 

The second mechanism which must be modeled for the utilization of transducers in control design is 
the magnetostriction due to changing magnetization. Once this model is obtained, it can be combined with 
the previous magnetization model to provide a characterization of strains output by the transducer in terms 
of currents input to the solenoid. For this investigation, we considered the quadratic model (3.1) as a first 
approximation to the relation between magnetization and magnetostriction. 

The performance of this model is indicated in Figure 5. At moderate drive levels, the strain data 
exhibits minimal hysteresis and is adequately characterized by the quadratic model. At the high drive 
levels illustrated in Figure 5b, the data exhibits significant hysteresis and saturates from a quadratic to 
nearly linear relationship as M approaches its maximum value. One component of this hysteresis is due 
to magnetostrictive hysteresis while other effects are due to mechanical hysteresis caused by the prestress 
mechanism. The performance of the model is much less accurate at high drive levels due to such unmodeled 



(a) 



(b) 


Figure 5. Experimental data ( ) and quadratic magnetostriction model dynamics (a) Low drive 

level, (b) High drive level. 
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dynamics. This loss in accuracy is also reflected in the change of the estimated saturation value = 
1221 x 10~ 6 in the low drive regime to X s = 1003 x 10 -6 at the high drive level. We note that at levels below 
that depicted in Figure 5a, the value A s = 1221 x 10~ 6 provides adequate model fits. 

The same tendencies are apparent when the magnetization and magnetostriction models are combined to 
provide a relationship between input currents and output strains. As illustrated in Figure 6a, the combined 
model is accurate at moderate drive levels and will be adequate for control design in this regime. Figure 6b 
illustrates that a high drive levels, however, the magnetostrictive model degenerates due to unmodeled 
dynamics and hysteresis. The extension of the magnetostrictive model to incorporate these effects at high 
drive levels is under current investigation. 




(a) (b) 


Figure 6. Experimental data (- ) and combined magnetization and magnetostriction model dynamics 

( ); (a) Low drive level, (b) High drive level. 


Minor Loop Model 

The modeling of minor asymmetric loops comprises the final component of this investigation. Accurate 
minor loop characterization is important for numerous applications including control design for transducers 
in unbiased and biased states. In a general unbiased state, it is crucial that the model be able to characterize 
both major and minor loop dynamics to attain the full range of dynamics specified by the control law. 
The characterization of minor loops in a biased state is important since it represents a common operating 
condition for transducers. 

For both cases, we employed the volume fraction (3.6) to attain closure in the minor loop magnetization 
model. The resulting model fit is illustrated in Figure 7a, b where 1.0 ksi data containing a major loop 
and two minor loops is considered. Both the major and minor loop dynamics are resolved by the model 
with the slight discrepancy in minor loop position due to differences in experimental and model major loop 
magnetizations for the values of H at which the turn points occur. We note that the model parameters 
for this case are those summarized in Table 2 and no parameter changes are necessary to accommodate the 
minor loops. As with the major loop case, turning points are dictated solely by the input magnetic field 
(or equivalently, the input current 7). Figure 7c and d illustrate the performance of the magnetization and 
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magnetostriction models in resolving major and minor loops in the 1.3 ksi data. The accuracy of the fit in 
Figure 7c reflects the accuracy of the underlying magnetization model while discrepancies in the major loop 
strain fit in Figure 7d are due to the previously mentioned unmodeled dynamics in the magnetostriction at 
high drive levels. At moderate levels, the minor loop model is sufficiently accurate for control applications. 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 0 0.5 1 1.5 2 2,5 3 

H {A/m) x10 4 H (A/m) * 10 4 


(c) (d) 

Figure 7. Experimental data ( ) and minor loop model dynamics ( ); (a) Two minor loops in magne- 

tization data with 1.0 ksi applied stress, (b) Magnified view of 1.0 ksi case, (c) Minor loop in magnetization 
data with 1.3 ksi applied stress, (d) Minor loop in strain data with 1.3 ksi applied stress. 

5. Concluding Remarks. An energy- based model for characterizing magnetization and output strains 
for magnetostrictive transducers is presented. The magnetization model, which is based upon the Jiles- 
Atherton mean field theory for ferromagnetic materials, provides a means of characterizing the magnetic 
hysteresis inherent to the transducer. Through enforcement of closure conditions, nested asymmetric minor 
loops as well as symmetric major loops are resolved by the characterization. This magnetization model is 
currently constructed for a transducer with quasi-static input and fixed operating temperature. Within this 
regime, the model provides the capability for characterizing variable input levels to the solenoid and differing 
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applied stresses to the Terfenol rod. The good agreement of this theory with experimental data illustrates 
the flexibility of the model under a variety of operating conditions. 

A quadratic model based upon the geometry of moment rotations was employed to quantify the mag- 
netostriction and strains generated by the transducer. As illustrated through comparison with experimental 
data, this characterization was adequate at moderate drive levels but degenerated at high drive levels due to 
unmodeled nonlinearities and hysteresis. Certain aspects of the magnetostriction hysteresis can be included 
through the energy model of [9] but adequate quantification of the full relation has not been attained and is 
under current investigation. 

At moderate input levels, the combination of the magnetization and magnetostriction models provide 
and accurate characterization of output strains in terms of input currents to the solenoid. For quasi-static 
applications in which temperature can be regulated, the model is sufficiently accurate for control design. The 
robustness of the model with regard to operating conditions and the small number of required parameters 
(six) enhance its suitability for such applications. 
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